Current-induced Vortex Motion by Spin- Transfer Torque 
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We investigate the dynamics of a magnetic vortex driven by spin-transfer torque due to spin 
current in the adiabatic case. The vortex core represented by collective coordinate experiences a 
transverse force proportional to the product of spin current and gyrovector, which can be interpreted 
as the geometric force determined by topological charges. We show that this force is just a reaction 
force of Lorentz-type force from the spin current of conduction electrons. Based on our analyses, 
we propose analytically and numerically a possible experiment to check the vortex displacement by 
spin current in the case of single magnetic nanodot. 
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Manipulation of nanoscale magnetization by electric 
current is one of the most attractive subjects in both 
basic physics and technological applications. After the 
theoretical prediction [1, 2], it has been widely recog- 
nized that spin-polarized current (spin current) plays a 
crucial role in magnetization dynamics. The spin cur- 
rent exerts a torque on localized spins by transferring 
spin angular momenta of electrons through the exchange 
interaction between conduction electrons and localized 
spins, which is thus called the spin- transfer torque. The 
key understanding of the effect is that the spin current 
favors magnetic configurations with spatial gradient, or 
more precisely, with finite Berry-phase curvature along 
the current. Such spatial gradient by spin current indeed 
gives rise to the motion of domain wall [3, 4], spin- wave 
instability in a uniform ferromagnet [5, 6, 7] and domain 
nucleation [8]. 

Recent experiments [9, 10] and numerical simulation 
[11] for current-induced domain wall motion have shown 
that there is a vortex-like configuration in magnetic 
nanowire. Also, magnetic vortices in nanodots have 
drawn much attention since the MFM observation of a 
vortex core [12]. However, an effective force on the vor- 
tex and its dynamics by spin-transfer torque due to the 
spin current have not been clarified. 

In this Letter, we present a microscopic theory of vor- 
tex dynamics in the presence of spin current by using the 



collective coordinate method. In the adiabatic approxi- 
mation, we derive an effective force exerted on the vortex 
core due to spin-transfer torque. It is shown that the vor- 
tex core experiences a transeverse force, which compen- 
sates the Magnus-type force derived from the so-called 
Berry's phase term. This specific force is topologically 
invarient, which is characterized by topological charges. 
Unlike the case of domain wall [4] , we show that there is 
no threshold current to induce the vortex motion in the 
absence of an external pinning. It is of great interest to 
the vortex-based devices for application. To verify the ex- 
isitance of current-induced transverse force on the vortex, 
we propose a possible experiment for the current-induced 
vortex displacement in a single magnetic nanodot. 

We consider the Lagrangian of the localized spins in 
the continuum approximation. The spins are assumed to 
have easy plane taken to be x — y plane, and are described 
by the Lagrangian 
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where S{x,t) = Sn{x,t) represents the localized spin 
vector with unit vector n = sin 9 cos 4> -t- sin 6 sin + 
coaO Ez and the magnitude of spin S; ei(i = x,y,z) 
are unit vectors of Cartesian frame. The J and K\ are 
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respectively the exchange and the hard-axis anisotropy 
constants and a is the lattice constant. The first term of 
the right hand side (r.h.s) in Eq.(l) is the so-called Berry- 
phase term, which in general determines the dynamical 
property of the localized spins. 

Let us denote the spin configuration of a vortex cen- 
tered at the origin by a vector field nv(a;) with imit mod- 
ulus. As a vortex profile, we take an out-of-plane vortex; 
ny{x — > 0) = pcz, where p = ±1 is the polarization, 
which refers to the spin direction of the vortex core center 
andn(|a;| ^ Sy) = cos{q(p+CTT/2)ex+sm{q(p+CTT/2)ey, 
where Sy is the vortex core radius, = tan~^{y/x), 
q = ±1,±2,..., is the vorticity, which describes the 
number of windings of the spin vector projected on the 
order-parameter space and C = ±1 is the chirality, 
which refers to the counterclockwise (C = 1) or the 
clockwise (C = —1) rotational direction of the spin in 
the plane. We here introduce a collective coordinate 
X{t) = X{t)e.j; + Y{t)ey, which represents vortex core 
center, and assume that a moving vortex can be written 
as n{x,t) = ny{x — X{t)) at least as a first approxima- 
tion, that is, ignoring the spin wave excitation. Substi- 
tuting this into Eq.(l), we obtain the Lagrangian for the 
collective coordinate as 



is the energy splitting, and cr is a Pauli-matrix vector. 
Here we perform a local gauge transformation [4, 13] 
in electron spin space so that the quantization axis is 
parallel to the localized spins n{x, t) at each point of 
space and time; c{x,t) — U{x,t)a{x,t), where a{x,t) 
is the two-component electron operator in the rotated 
frame, and U{x, t) = m{x, t) ■ cr is an SU(2) matrix with 
m = sin(0/2)cos0 Sx + sin(0/2) sin (/) By -\- cos(^/2) e^- 
The Lagrangian is now given by 



ih (do + iAo) 



A / d^a;aV^a, (6) 



Ly = ^G-{X xX)- U{X). 
Here G is the gyrovector defined by 



(3) 



where = A,y ■ a = -iW 8,^11 {v = 0, x, y, z) is SU(2) 
gauge field determined by the time and spatial deriva- 
tive of the localized spins. For slowly varying magnetic 
configurations, the electron spins can mostly follow adia- 
batically the localized spins. This is justified for the con- 
dition [Is./ e-p)k-p\ ^ 1, where ep and fcp are, respectively, 
the Fermi energy and Fermi wave number of conduction 
electrons, and A is the characteristic length scale of the 
spin texture of the localized spins. In this adiabatic ap- 
proximation, taking the expectation value of Lei for the 
current-carrying nonequilibrium state, we can obtain the 
following interaction Hamiltonian of the first-order con- 
tribution to the localized spins [8]; 



/dPx hS f 
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with L being the thickness of the system, and U{X) is 

a potential energy of a vortex evaluated from the Hamil- 
tonian Hg. The gyrovector G is topologically invariant 
corresponding to polarization p and vorticity q and the 
number of corvering the mapping space D. In the case 
of nv{x), this mapping number is 1/2 in unit of surface 
area iw. 

The first term of the r.h.s in Eq.(3), which comes from 
the Berry-phase term, represents that X and Y are essen- 
tially canonically conjugate each other. This term pro- 
vides a transcverse force —GxX on the moving vortex, 
the so-called Magnus force, perpendicular both to the gy- 
rovector and to the vortex velocity, whose term has been 
derived and discussed by many workers [15, 16, 17, 18, 19] 
in various systems. 

Let us investigate the force acting on the vortex by 
spin current of conduction electrons. The Lagrangian of 
the electrons are given by 



where js is the spin current density, which is written by 

using the distribution function fka = {o-l^a-'^kcr) iu the 
rotated frame specifying the current-carrying nonequi- 
librium state as jg = y ^ a^fk^.As seen from Eq. 
(7) , the spin current favors a finite Berry-phase curvature 
along the current. Indeed, the Hamiltonian Hst leads to 
the spin- transfer torque in the case of domain wall [8] . 

From Eq.(7), we can derive the effective Hamiltonian 
represented by the collective coordinate for the vortex 
core. Taking the variation 5n = —djny{x — X)6Xj, 
where repeated roman indices imply sum over the in- 
plane spatial direction j = x, y, we obtain 

6HsT = ^2TTLpq{v,,xSY - v,,y6X), (8) 



d^xc\x,t)Un— + —V^ic{x,t)-H,d, 



(5) 



where c(c^) is anihilation (creation) operator of conduc- 
tion electrons. The last term Hgd represents the exchange 

interaction between localized spins and conduction elec- 
trons given by ffgd = —A J d^xn ■ (c^crc)x, where 2 A 



where Vs = represents the drift velocity of electron 

spins. By integrating Eq.(8), we obtain 

HsT = G-{vsXX). (9) 

Thus, a force acting on the vortex core is given by 

(10) 



dX 

This current-induced transverse force has been previ- 
ously derived by Berger [14] based on a phenomenological 
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treatment in the case of Bloch line. Here we have derived 
this force from microscopic theory. It is noted that this 
force does not depend on the chirahty C = ±1 of the vor- 
tex in contrast to a force produced by a magnetic field. 
Since it is hard to control the chirality, this fact would 
great advantage in application. A microscopic deriva- 
tion of a general relation between force and torque in the 
Landau-Lifshitz-Gilbert (LLG) equation is presented in 
Ref. [20]. 

Before proceeding, we here briefly remark that this 
force can be interpreted as a reaction force of a Lorentz- 
type force from the spin current of conduction elec- 
trons. The interaction Hamiltonian, Eq.(7), can be 
rewritten by using the the U{1) gauge field, A{x) = 
■^V(t){x){cos9{x) — 1), interacting with the spin current 
density. Thus the magnetic field B{x) is given by 



B = {da:Ay-dyA^,)e:, 



2e 



n-{d^ny.dyn)e^, (11) 



which corresponds to the so-called topological field [21]. 
It is noted that the SU(2) field intensity, F^i, = d^Ai, — 
duA^ — i[A^j^, Ay]^ vanishes, since the original Lagrangian 
L[^[ in Eq.(5) docs not include the local gauge field. The 
finite magnetic field is a consequence of the projection 
from SU(2) to the U(l) by taking the adiabatic approx- 
imation. Under this magnetic field, the spin current of 
the conduction electrons may experience the following 
Lorentz-type force as 



d X js X B = G X Vs 



^ST, 



(12) 



which can be interpreted as the reaction force acting on 
the vortex. 

Let us derive the equation of motion for the collective 
coordinate of vortex in the presence of spin current based 
on the Euler-Lagrange equation; 



d dL dL _ dW 



(13) 



where L — Ly — Hst is the total Lagrangian, W is the 
so-called dissipation function given by 



hS f d^x .9, , 01 
W^a— n\x,t)^-DX\ 



(14) 



with a being the Gilbert damping constant. The con- 
stant 

Pi Q r 

D = ^L dxdy {{di9f + sin^ e{di(t)f] , (15) 
a Jd 

generally includes a factor ln(i?v/'5v), where i?v is the 
system size. Here we assume that system has rotational 
invariance along the z axis. The concrete expression in 
Eq. (13) is given by 



Gx{v^- X) 



dU{X) 
dX 



aDX. 



(16) 



This is the equation of motion for the vortex dynamics 
in the presence of spin current. If the r.h.s in Eq.(16) is 
absent, we obtain X — Vg, where the vortex core moves 
along the spin current perpendicular to the transverse 
force -FsT- This situation can be seen from the current- 
induced domain wall motion [4, 8]. On the other hand, 
the damping term —aDX acts as deviating from the 
orbital direction of the moving vortex along the current. 
Importantly, there is no intrinsic pinning in the dynamics 
of vortex unlike the case of domain wall. This lead to a 
vanishing threshold current for the vortex motion in the 
absence of an external pinning. This is because, in the 
translationally-invariant system, there is no pinning on X 
and Y. This is in contrast with the case of domain wall, 
where, 0o is pinned by the hard-axis magnetic anisotropy 
even in the translationally-invariant system [4]. Thus 
vortex-based devices would have great advantage in low- 
current opperation. 

To verify the exisitence of current-induced force on the 
vortex, we propose the vortex displacement by spin cur- 
rent in single magnetic nanodot [22], where an out-of- 
plane vortex with vorticity q = 1 is stabilized. We assume 
the electric current is uniform in the nanodot, and fiowing 
in the possitive x-direction, that is, Vg — -^js — 'Vs^x- 
We assume the full spin polarization of the current, 
P — 1, for simplicity (Fig. 1). The potential energy 
U{X) is modeled by a harmonic one U{X) = kX^/2, 
where k is a force constant. In Ref. [23], k, is evaluated 
in detail, which depends on the aspect ratio g = L/R, 
where R is the dot radius. From Eq.(16), the equation of 
motion is given by 



(1 + ia)Z — —iLoZ + Vs 



(17) 



where Z ^ X + iY, a = aD/G, and uj = k/G. For an 
initial condition Z[0) = 0, the solution is given by 



Z{t) = i— <{ exp ( - /^^ , 



(18) 




FIG. 1: Schematic illustration of an experimental set up for 
current-induced vortex displacement. The topological charges 
are chosen to be p = 1, g = 1, and C = 1 in the above vortex. 
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Thus the vortex center exhibits a spiral motion, whose 
rotational direction depends on the sign of the core po- 
larization p = ±1 (Fig. 2(a)). The final displacement of 
the vortex core is perpendicular to the current direction 
and given by 5Y = —Gvs/k, which also depends on p, 
where the transverse force —GxVg balances the restor- 
ing force —kX. This result is consistent with the recent 
experimental one in Ref. [10], where a distorted vortex- 
wall was shifted to the diagonal direction and stopped 
moving. 

For comparison, the current-induced dynamics of a 
vortex was calculated by micromagnetic numerical sim- 
ulations based on the LLG equation with spin current 
terms; 

dM „ a dM 

— = --foMxH,ff + —Mx — -v,-VM, 19 
ot Als at 

where M{x,t) is the magnetization vector, 70 is the gy- 
romagnetic constant, Hcs is the effective magnetic field 
including the exchange and demagnetizing field and Ms is 
the sarturation magnetization. The last term represents 
the spin-transfer torque [5, 6, 7, 8]. The sample is di- 
vided into identical cells, in each of which magnetization 
is assumed to be constant. The dimension of the cells is 
4 X 4 X /i nm^ with dot thickness h = 10, 20, 30 nm. 
The dot radius is taken to be i? = 500 nm. The com- 
putational material parameters are typical for permal- 
loy; Ms = 1.0 T {Ms/fio = 8.0 X 10^ A/m), the ex- 
change stiffness constant A — 1.0 x 10^^^ J/m and 
70 = 1.8 X lO^m/A • s. We take a = 0.02. The pro- 
graming code is based on those of Rcfs. [11, 24]. 

Figure 2(b) shows the time evolution of volume- 
averaged magnetization, which exhibits the spiral motion 
of vortex core. It is noted that the rotational direction 
is opposite to the case of analytical result because of the 
replacemet TiS/a^ —Ms/70. Figure 3 shows that vor- 
tex displacement as a function of the spin current density 
in various aspect ratios g — L/ R. The numerical results 
are in good agreement with the analytical ones for small 
vortex displacement. It is found that the smaller g is 
more advantageous to the vortex displacement. 




In conclusion, we have clarified the transverse force on 
the vortex and its dynamics by spin-transfer torque due 
to the spin current in the adiabatic regime. We have pro- 
posed analytically and numerically a possible experiment 
for the vortex displacement by spin current in the case 
of single magnetic nanodot. Finally, this vortex displace- 
ment may affect a magnetoregistance, which is probable 
to be detected, for example, by using planar Hall effect 
[25]. 

The authors would like to thank T. Ishida and T. 
Kimura for valuable discussion. 
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